Dissipation and detection of polaritons in ultrastrong coupling regime 
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We have investigated theoretically a dissipative polariton system in the ultrastrong light-matter 
coupling regime without using the rotating-wave approximation on system-reservoir coupling. Pho- 
tons in a cavity and excitations in matter respectively couple two large ensembles of harmonic 
oscillators (photonic and excitonic reservoirs). Inheriting the quantum statistics of polaritons in the 
ultrastrong coupling regime, in the ground state of the whole system, the two reservoirs are not in 
the vacuum states but they are squeezed and correlated. We suppose this non-vacuum reservoir 
state in the master equation and in the input-output formalism with Langevin equations. Both two 
approaches consistently guarantee the decay of polariton system to its ground state, and no photon 
detection is also obtained when the polariton system is in the ground state. 

PACS numbers: 42.50.Pq, 03.65.Yz, 42.50.Lc 



I. INTRODUCTION 

Light-matter ultrastrong coupling [1-28] means that 
the coupling strength g is comparable to or larger than 
the transition frequency w x of excitations in matters 
(g > lu x ), and it shows a variety of peculiar proper- 
ties, such as virtual photons in the ground state [1, 5], 
squeezed eigen states [1, 19], nearly degenerated ground 
states [11, 17, 18], quantum phase transitions [12, 29], 
and so on. The ultrastrong coupling has been realized 
experimentally by intersubband transitions in semicon- 
ductor quantum wells [6-9, 26, 27], artificial atoms in su- 
perconducting circuits [14-16], and cyclotron transition 
in two-dimensional electron gas [25]. In most cases, pho- 
tons with THz or microwave frequency are confined in a 
cavity, and the cavity photons are coupled with external 
photonic field (outside the cavity) by unignorable dissi- 
pation rate (but small compared to the light-matter cou- 
pling). In other words, when we neglect the coupling with 
matters, the real eigen modes of photons are represented 
as coupled fields of cavity mode and external fields [30]. 
As we will discuss in the present paper by the Fano-type 
diagonalization technique [31-33], when the cavity mode 
is squeezed in the ground state due to the ultrastrong 
light-matter coupling, the external photonic field is also 
squeezed in the ground state of the whole system. How- 
ever, we cannot observe the squeezing nor energy flow by 
photon detectors if the system is in the ground state. In 
this paper, we have developed two frameworks, the mas- 
ter equation and Langevin equations with input-output 
relation. They are derived by supposing the squeezed ex- 
ternal fields and show no photon detection consistently. 

In the standard theory of quantum optics [2, 24, 34- 
42], in order to introduce a dissipation of cavity mode 
or of excitations in matters, we consider coupling with 
an ensemble of harmonic oscillators with continuous fre- 
quencies, and the oscillators are supposed to be in the 
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vacuum state. The dissipation of focusing system has 
been successfully described by such a treatment in both 
master equation and input-output formalism, at least 
in the weak and (normally) strong light-matter coupling 
regimes (g < r and r < g <C cj x , respectively, for dissipa- 
tion rate T). As pointed out in some papers [34, 39-43], 
the master equation should be derived by considering the 
eigen states of the focusing system, and the rotating-wave 
approximation (RWA) should be performed carefully on 
the system-reservoir coupling even if the system-reservoir 
coupling is weak compared to the light-matter coupling 
(in the strong light-matter coupling regime). In the ul- 
trastrong coupling regime, such treatment has been per- 
formed by Beaudoin, Gambetta, and Blais [24]. The dis- 
sipation of the ultrastrong coupling systems can be suc- 
cessfully described under the RWA on system-reservoir 
coupling (both pre-trace and post-trace RWAs are used 
in terms of Ref . [41] ) by considering the eigen states of the 
cavity system (there are squeezed virtual photons in the 
ground state). Furthermore, the photon detection has 
also been discussed in Ref. [28], and virtual photons in 
the cavity are not counted by normal- and time-ordering 
the operators of cavity system. However, the RWA on 
the system-reservoir coupling destroys the information 
on how the cavity system couples with the reservoirs. In 
other words, the reservoirs become a black box, and we 
effectively suppose simplified reservoirs coupling with the 
eigen states. Although such a treatment is appropriate 
and simple for discussing the dynamics of the focusing 
system, it becomes difficult to discuss the statistics of 
output photons emitted from the cavity. There are ac- 
tually virtual photons not only in the cavity but also in 
the photonic reservoir even if the whole system is in the 
ground state, which is naturally derived by the Fano- 
type diagonalization, although the virtual photons are 
not counted by detectors. Of course, when the cavity 
system is excited, we can detect photons emitted from 
the cavity. Whereas antibunching of emission can sur- 
vive even under the RWA on system-reservoir coupling 
[28], the quantum fluctuation (squeezing) of the emis- 
sion is easily diminished in such treatment, because the 
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interference between reservoir free field and cavity contri- 
bution is important for squeezing [35, 36, 38]. Therefore, 
in order to fully discuss the quantum statistics of emis- 
sion in the ultrastrong light-matter coupling regime, we 
have to develop a comprehensive framework describing 
the dissipation and emission without using the RWA on 
light-matter coupling nor on system-reservoir one. 

In the weak and normally strong light-matter coupling 
regimes, the ground state of the whole system is approxi- 
mately represented by the vacuum states of photonic and 
excitonic reservoirs (exactly the vacuum states under the 
RWA on light-matter coupling). In both master equa- 
tion and input-output formalism, the vacuum reservoirs 
are usually considered for describing the dissipation, and 
no photon detection is naturally obtained in the ground 
state. Then, these three approaches (analysis of ground 
state, master equation, and input-output formalism) are 
consistent in the standard dissipation theory in the weak 
and strong light-matter coupling regimes. However, if 
we simply suppose the vacuum reservoirs, in the ultra- 
strong light-matter coupling regime, the master equation 
and input-output formalism give different results. This 
is because the photonic and excitonic reservoirs are not 
in the vacuum states in the ground state of the whole 
system, but they are actually squeezed and correlated. 
We have to suppose this non-vacuum reservoir state in 
order to remove the discrepancy between the results of 
master equation and of input-output formalism. 

In the present paper, we discuss polariton system con- 
sisting of two bosonic modes, photons in a cavity and 
excitations in matter, each of which couples with an en- 
semble of harmonic oscillators (photonic and excitonic 
reservoirs). Diagonalizing the whole system by the Fano- 
type technique [31-33], we find the squeezed and corre- 
lated reservoirs in the ground state. By supposing this 
reservoir state, the master equation certainly guarantees 
the decay of the polariton system to its original ground 
state in the closed case. We also check that, when the 
polariton system is in the ground state, the virtual pho- 
tons in the photonic reservoir are not counted by normal- 
and time-ordering the operators in polariton base. In the 
input-output formalism, we also obtain no photon detec- 
tion by supposing the squeezed and correlated reservoirs 
and by normal- and time-ordering the operators. Then, 
we achieve the consistency of the three approaches (di- 
agonalization, master equation, and input-output formal- 
ism) even in the ultrastrong light-matter coupling regime. 

This paper is organized as follows. The Hamiltonian 
is shown in Sec. II, and basic features of the ultrastrong 
coupling regime is also discussed. The Fano-type diago- 
nalization of the photonic part is performed in App. A. 
In Sec. Ill, we diagonalize the whole Hamiltonian and 
show that the reservoirs are squeezed and correlated even 
in the ground state of the whole system. The master- 
equation approach is discussed in Sec. IV, where we 
demonstrate the decay of polariton system to its original 
ground state. Correlation functions of reservoir fields are 
calculated in App. B, and the detailed calculation of mas- 



ter equation and photon detection is shown and App. C. 
The input-output approach is discussed in Sec. V, and the 
detailed calculation of photon detection in this approach 
is shown in App. D. Finally, we discuss the comparison 
with previous theories in Sec. VI, and the summary is in 
Sec. VII. 



II. HAMILTONIAN 

The Hamiltonian describing cavity photons and exci- 
tations in matter is written as 

H s = hoj c a^a + hu x tfb + ihg(a + a^)(b-tf) + HD(a + a^) 2 . 

Here, a and b are annihilation operators of the photon 
and excitation, respectively, satisfying the bosonic com- 
mutation relation [a, a'] = [b , b'] = 1 and [a, a] = [b, b] = 
[a,b] = [a, 6 T ] =0. lo c and w x are their eigen frequencies, 
and g is the coupling strength. The ultrastrong coupling 
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The last term is the so-called diamagnetic 
term naturally derived in the minimal coupling scheme 
[1], and the coefficient is normally D > g 2 /uj x , by which 
we cannot expect the quantum phase transition [12, 29]. 
If the polariton system is isolated from the environment, 
as discussed in Ref. [1], this Hamiltonian can be diago- 
nalized as 



Hs = h^jP^jPj + const. 

j=L,U 



(2) 



Here, ph and pu are annihilation operators of lower and 
upper polaritons, respectively. They are represented as 
a combination of annihilation and creation operators of 
photon and excitation: 



Pj = Wjti + Xjb + yja) + zj 



(3) 



These coefficients and cigen frequencies uij are deter- 
mined by solving 
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V -ig ig -cj x 

(4) 

From this eigen value problem, we get four eigen values 
{u>L,LJu, —ujl, — <-<Ju}i whose eigen vectors correspond to 
{pl,Pu,PliPu}' respectively. The coefficients are nor- 
malized for satisfying [Pj,pt] = Sj t k for j, k = L,U and 
we also get [pj,Pk] = 0. Inversely, the photon and excita- 
tion operators are represented by the polariton operators 
as 
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The photonic and cxcitonic reservoirs are individually 
represented as ensembles of harmonic oscillators as 



(G) 



where a m and (3 m are annihilation operators of oscilla- 
tors in photonic and excitonic reservoirs, respectively, 
and i?^ x is the oscillating frequency. The ensembles 
show nearly continuous spectra. These operators satisfy 
[d m ,d4] = \j3 m ,Pn\ = Sm,n and [«m,a„] = \fi m ,$n] = 
[d m ,/3„] = [a m ,/3t] = 0. The system-reservoir coupling 
is represented as 

F s -R = - a^F c ) + ih(Pfb - PP X ). (7) 

Here, F c and F x are photonic and excitonic reservoir 
fields, respectively, and they are expressed by the an- 
nihilation operators a m and (3 m and coupling strengths 
and 7 m as 



n i 



(8b) 



It is worth noting that Eq. (7) is not the result of 
RWA, but this expression is naturally derived consid- 
ering the transmission and reflection of particles be- 
tween inside and outside of the cavity [30]. In other 
words, concerning the coupling between cavity photons 
and external photonic field, they are coupled through 
the electric field and also through the magnetic field. By 
summing these two interaction (ifi,/2)(a — a)){F^. + F c ) 
and (ih/2)(a + at)(Ft — F c ), we can derive the first 
term in Eq. (7). For simplicity, we also suppose the 
similar situation concerning the coupling with excitonic 
reservoir. Whereas the Hermitian expressions has been 
supposed for describing the dissipation in some works 
[24, 28, 39, 40], Eq. (7) can be considered as the stan- 
dard expression, because there is no ambiguity whether 
the system-reservoir coupling is electric or magnetic. Of 
course, when we suppose specific systems, the expression 
of system-reservoir coupling is automatically determined. 
In terms of the polariton operators, the system-reservoir 
coupling is rewritten as 

Hs-r = HPIpl - p[P L ) + \K{fIpu - p\jFu), (9) 
where 



F S = w 3 F c 



Ft 



%Ft 



(10) 



is the external field that couples with the lower (j = L) 
or upper (j = U) polariton inside the cavity. 

As discussed in Ref. [1], when we consider only the 
Hs system isolated from the reservoirs, there are virtual 
photons and virtual excitations even in the ground state. 
This is because the ground state should satisfy pj |g) = 0, 



and then the cavity mode is represented by a squeezed 
vacuum state in the ground state as 



(aa) g = - u 

j=L,U 



(11a) 
(lib) 



where (. . .) g means an expectation value in the ground 
state |g). The excitations in matter are also expressed 
as a squeezed vacuum state, and the photons and exci- 
tations are correlated in the ground state as discussed in 
Ref. [1]. If the system-reservoir coupling is weak enough 
compared to the light-matter coupling g, the squeezing 
and correlation of cavity photon and excitation should 
be maintained. 



III. DIAGONALIZATION OF WHOLE SYSTEM 

First of all, we diagonalize the whole system H = Hs + 
Hs-r + Hr by using the Fano-type technique [31-33]. As 
shown in App. A, the photonic part consisting of cavity 
mode and photonic reservoir is diagonalized as 

H ph = huj c a^a + £ htt^atlnOtm + ih(F%a - a^Fc) 

m 

(12a) 

) 

dw ftijA^ (u)A(u>) + const. (12b) 



Here, the partially diagonalized operator A(tu) is defined 
in Eq. (A3) and satisfies 



A(u;),H ph = HuiA(uj) 



and 



(13) 



(14) 



A(u),A\u>') ^S(lu-lo'). 
Similarly, the excitonic part is diagonalized as 

H cx = hu x b r b + £ hQplPj + ih(Flb - PP X ) (15a) 



du! hujB ] (uj)B(uj) + const. 



(15b) 



The operator B(lu) is represented in Eq. (A12). Using 
these partially diagonalized operators, the whole Hamil- 
tonian is represented as 



H 



du [tvuti (w)A{u) + HljB^(lj)B(ui)} 



+ \hg{a + a f )(S - S f ) + HD(a + a f ) 2 + const. (16) 



The light-matter coupling and diamagnetic terms arc also 
expressed in terms of A(lu) and B(uj) by using Eqs. (All) 
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and (A13). For diagonalizing the whole Hamiltonian H, Then, H can be diagonalized as 
we suppose a new operator 



H 



dto fiwF(w)P(w) + const. 



(20) 



P(co) = / du>' [W{lo,lo')A(lo') + X{lo,lo')B(oj') 
Jo 

+ Y{lo, tu')A\uj r ) + Z{lo, lo')B\lo')]. (17) 
The coefficient functions are determined for satisfying 

[P(oj),H]=tlu)P(u}) (18) 

and 

From Eq. (18), the coefficient functions are determined 
[P(ui), P^(lo)] = S(lo — lo'). (19) by the eigen value problem 



/lo' + 2D-lo -ig -2D -ig \ / u c {w')W{w,w') 

ig lo' — lo —ig I u x (lo')X(lo,lo') 

2D -ig -lo'-2D-lo -ig u c (u)')*Y(oj,u)') 

\ —ig ig —lo'—lo/ \u x {lo')* Z{lo,lo')/ 



(21) 



where coefficients u c ^ x {lo) is represented in Eq. (A9). This eigen value problem is equivalent to Eq. (5) by replacing 
lo c and lo x with lo' , and the eigen frequencies lol,u{u') must be equal to lo. Then, operator P{lo) is represented as 



(22) 



Here, lo'j is the frequency satisfying lo = u)j(u)'j), and where coefficient v c (to', lo) is expressed in Eq. (A6). The 
Wj(io), Xj(co), Dj(uS), and Zj(oj) are the coefficients when phase-independent correlation function is written as 



we solve Eq. (5) by replacing to c and lo x with lo. They 
are normalized for satisfying Eq. (19). Inversely, we can 
rewrite A(co) and B(lo) as 

A(lo) = Uc (lo) [«>iM*-P(wjH) - ^H^KM) 

j=L,U 

(23a 

j=L,U 

(23b) 

Then, the original photon, excitation, and reservoir op- 
erators are also expressed in terms of P{lo) by using 
Eqs. (All) and (A13). 

The ground state |g) of the whole system H is deter- 
mined for satisfying P(w)|g) = for < lo < oo. Then, 
when we apply the photonic free field 

Fj°)(r) = jH R tp(p) c -iH R t ( 24 ) 
onto the ground state, it is represented as 
# c (°)(r)|g) = - / dui/ da/ e-' 1UT K{io)v c {u',io)*u c {uj') 

(25) 



o 



' k(lo)* k(lo') 



dto I dco' 

IQ JO JQ 

xv c (lo",lo)v c (lo",lo')*\u c (lo")\ 2 £ \y 3 (io")\ 2 . (26) 

j=L,U 



Here, the coefficients u c {lo) and v c {lo,lo') are singular at 
lo = co c as seen in Eqs. (A9) and (A6). Since the coeffi- 
cients are normalized as 



dto \u c (lo)\ 2 = 1, 



(27) 



dLo v c {lo,lo')v c {lo,lo") =S(u)' -lo"), (28) 



if the dissipation of photons is weak enough compared 
to the characteristic frequencies of H$ system (g, lo c , 
and w x ), the correlation function is approximately repre- 
sented in the ground state as 

(g|Fi )tFi°)(r)|g) 

dLo c~^ t \k(lo)\ 2 l%K>| 2 = G c {r)(tfa) s . 

j=L,U 

(29) 
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In the same manner, the phase-sensitive correlation is 
expressed in the ground state as 

<g|Fj°>jK°>(r)|g) 

rOQ pOO pQC 

= - dcj duj' duj" c-' 1u) ' t k{uj)k{uj') 
Jo Jo Jo 

xv c (Lo",aj)*v c (u J ",Lu'yu c (u J ") 2 WjWYViW) 

j=l,u 

do; e~ iuT \K(uj)\ 2 Y, w (oo c yv 3 {io c ) 
= G c {T)(aa) g . (30) 



Therefore, even in the ground state, the photonic reser- 
voir field F c is also squeezed by the same order as the 
cavity mode. In the same manner, the correlation of the 
two reservoir fields are also the same as internal ones 
if the dissipation is weak enough. In other words, the 
internal modes a and b are balanced with F c and F x , 
respectively, in the ground state of the whole system. 

When we initially suppose a squeezed cavity mode and 
a big reservoir in the vacuum state, of course the reservoir 
does not become equally squeezed but it almost remains 
in the vacuum state after switching on the coupling be- 
tween them. This is a non-equilibrium problem. How- 
ever, when we consider the equilibrium of the cavity mode 
and the reservoir, if the ground state of the cavity mode 
is squeezed, the reservoir field is also squeezed in the 
ground state of the whole system. Therefore, when we 
consider the dissipation of the ultrastrong light-matter 
coupling system to its original ground state, we should 
suppose the squeezed and correlated reservoirs as will be 
discussed in the following two sections and also in Sec. VI. 



IV. MASTER-EQUATION APPROACH 

In this section, we derive a master equation for de- 
scribing the dissipation of polariton system Hs by con- 
sidering the coupling with reservoirs. Obeying the stan- 
dard derivation of master equations [24, 37, 39-42], from 
the expression (9) of system-reservoir coupling -ffs-R; the 
master equation for reduced density operator p(t) de- 
scribing Hg system is derived under the Born approxi- 
mation as 



where £ — Co + £diss and 



Co\p] = T^[Hs,p{t)] 



,P 



Pj,Vfk\p] 



j,k=L.U 

+ [cf k [p],P 3 ] + [p},Cfk[p]\ + [p],Bj k \p] 
Bfk[p},p]\ + [pj,^ k [p]\ + [J^ k \p],Pi]} 



(33) 



J to 

(34a) 

T>l[p] = fdt' Mt-t')W)pl](Fi \t')Ff\t))- m 



(34b) 



[p] = ?7 s (t-i')[pi/5(t')](^ (0)t (i)^ 0) (i'))i 

Jto 



(34c) 



cf k [ P ] = I m Mt-*)W)pM 0) \*)F?\t))i 



to 



(34d) 

in 

(34e) 



Bf k [p}= / dt't> s (t-tO[p(i')^](4 (0) (*')^ 0) (t)>i 

Jto 



(34f) 



A) k [p] = I dt'U s (t-t')\p k p(t%Fj 0) ^t)Fl 0) \t% 



to 



(34g) 



Af k [p]= / dt' U s (t-t')W)M(Ft 0) Ht')Fl 0) \t)) 



(34h) 

Here, to — > — oo is the switch-on time of system-reservoir 
coupling. Us(t)[0] is the propagator in Hs system for 
arbitrary operator O as 



U S (t)[0) = c ~iH s r/hQ c iH s r/h^ 



(35) 



and F^jj is the reservoir field in the interaction picture 
(free field): 



F^ht) = e i6Rt F LtU e- iARt . 



(36) 



(31) These free fields are in the polariton base, and it is rep- 
resented by the ones F^x in the excitation-photon base 
as in Eq. (10). They satisfy (/x, v = c, x) 

[Fj l \t),FW(t')} = 6 lltV G ll (t-t'), (37a) 

(32) [PM(t),FP(t')} = 0. (37b) 
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Time t / (2ti/coc) 



FIG. 1: Numbers of lower and upper polaritons are plotted 
as a function of time t. The initial state is given as p(0) = 
|g)(g|, and the reservoirs are supposed to be in the vacuum 
state. The time-development is calculated by master equation 
(31), correlation (39), and memory kernel (40). Parameters: 
w x = lo c , g — lu c , D = g 2 /uj x , r c = r x = lCT 2 aj c , and 

^cut-off /-icut-off' i riS 
ll c = H x =10 OJc. 

where G c (t) and G x (t) are memory kernels due to the 
coupling with reservoirs and arc expressed as 

G c (T)=J2\K m \ 2 e- in ™ T , (38a) 
m 

G x (r)=J2hm\ 2 e- ia ^. (38b) 

m 

Prom the master equation (31), the dynamics in Hs sys- 
tem are determined by supposing correlation functions of 
the free 1 fields, which arc 1 considered as an input from the 
reservoirs to the H$ system. 

First of all, we suppose that the photonic and excitonic 
reservoirs are in the vacuum state and the correlation 
functions are given as (fi, v = c,x) 

(F^(t)F^(t')) in = S^G^t - f), (39a) 
(F^(t)F^(t% n = (^ 0) (t)F(°)(t'))in = 0. (39b) 

In Fig. 1, supposing the ground state p(0) = |g)(g| at the 
initial time t = 0, we plot the development of numbers 
of lower and upper polaritons calculated by the master 
equation (31) and the correlation (39). The memory ker- 
nels are simply supposed as 

G c (t) = J°° dQ §c" ifiT , (40a) 

^cut-off 

G x (r) = f X dtt ^e~ inr , (40b) 



where the cut-off frequency f2^ governs the memory 
time of the reservoirs as ~ 1/ J?^ t_off . In a similar way as 
in Ref. [5] , the density operator p{t) is moved outside the 
time integral in the master equation (31). This treatment 
is valid if the memory time 1/Q^:~° S is short enough 
compared to the specific oscillation periods (l/w c , l/ui x , 
and 1/g) of Hs system. 

As seen in Fig. 1, the polaritons are excited by the 
vacuum reservoirs (at zero temperature). The periods of 
the oscillation are approximately tt/lul.u {^l/^>c = 0.414 
and lju/oj c = 2.414), whereas they are slightly modi- 
fied by the Lamb shifts. After a long time compared 
to 1/T C = 1/T X = (100/2tt) x (2tt/u) c ), the numbers of 
polaritons reach to certain values, which depend on the 
system-reservoir coupling strengths F c x . 

The polariton system is excited by the vacuum reser- 
voirs, because it is excited when the virtual photons in 
the ground state escape to the reservoirs. In other words, 
the ground state of polariton system is modified by the 
coupling with vacuum reservoirs. This result can also be 
understood by Eq. (33). In order to guarantee the de- 
cay of the Hs system to its original ground state |g), the 
photonic and excitonic reservoirs should not be in the 
vacuum state in the excitation-photon base (in terms of 

Fc°]), but the free fields F^ in polariton base should be 
in the vacuum state. In Refs. [24, 28], owing to the RWA 
on system-reservoir coupling, the decay to the ground 
state |g) is guaranteed by simply considering the vac- 
uum reservoirs in the excitation-photon base. However, 
in the present paper, we do not use the RWA to maintain 
the information of quantum fluctuation of the reservoirs. 
Instead, we suppose that the reservoirs are in the vac- 
uum state in polariton base (squeezed and correlated in 
excitation-photon base). 

Let's derive the correlation of reservoir free fields F { °1 

that guarantees the decay to the ground state |g) of Hs 
system and is simultaneously appropriate to the anal- 
ysis of Fano-type diagonalization discussed in Sec. III. 
We assume that the Hs system is in the ground state 
as p = |g)(g|. Under this assumption, let's inversely con- 
sider how the reservoirs are modified by coupling with Hs 
system. As discussed in Ref. [35] and in App. B of this 
paper, we can derive the correlation of free fields Fj^\j 

(on output side) from the density operator p{t) of Hs 
system. The equations of motion (Langevin equations) 
of cavity photons and excitations are derived as 

|U(t) = 1[S, H s ](t) - J" df G c (t - t')a(t>) - F^(t), 

(41a) 

|&(t) = ^[b,H s ](t) J* df G x (t t')b(t') F(°\t). 

(41b) 

In the standard theory of quantum optics [36, 38], the 
memory kernels G c ^ x {t) arc approximately described by 
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the Dirac's delta function by elongating the frequency 
range of reservoirs to — oo and oo. Since the escaped 
photons do not reenter into a cavity, wc can consider 
that the photonic reservoir has a quite small coherence 
time, and this approximation seems valid in most cases. 
However, in the ultrastrong coupling regime, when we 
do not use the RWA on system-reservoir coupling, we 
must keep the reservoir frequencies positive [2], and the 
Langevin and master equations are written in the time 
nonlocal forms in general. Although in the case of time- 
local equations we usually use the standard input-output 
relation [36, 38], we calculate the correlation between the 
free fields and internal ones by the formalism of Ref. [35] . 
First, we define the propagator U(t)[- ■ ■ ] satisfying 



is negligible compared to Eqs. (46). Furthermore, the 

self-correlation of free fields F c ,x is obtained in a steady 
state (the ground state in the present case) as 

(F^(t)F^)=G;(t)(sI§ v ) s (48a) 

(^ 0) W^ 0)t >=G,(r)(^4)g (48b) 

p(0)\ = jG>(|T|)(a M S„) g forr>0 

G„(\t\)(s^) s forr<0 1 ' 



p(t + r) = U(r)[p(t)h 



U(t) = C[U(t)} 



(42) 



(43) 



and the quantum regression theorem [35-38, 44, 45] is 
written for r > as 

<Oi(i + r)6 2 (t)> = Trid&iTftdiftt)]}, ( 44a ) 
{di{t)6 2 {t + t)> = Tr{6 2 U(T)[p(t)6 1 ]}. (44b) 

As discussed in detail in App. B, by using this and 

p(0) 



where (• ■ -) g means an expectation value in the steady 
state (ground state) . These correlation certainly satisfies 
Eqs. (37), and is equivalent to the ones derived in Sec. III. 

In the sense of perturbation theory, they are the corre- 
lation on output side, i.e., the modification of the reser- 
voirs due to the coupling with Hs system. The free-field 
correlation appearing in the master equation (33) is the 
one on the input side (effect from reservoirs to Hg sys- 
tem). Here, under the equilibrium between H§ and reser- 
voirs, the correlation of Fc°x should be equivalent on both 
input and output side. Let's substitute Eqs. (48) to the 
master equation (33). The frcc-ficld correlation in the 
polariton base can be derived by using Eq. (10). From 
Eqs. (48), we can easily get for j, k = L,U 

<iH 0) (r > 0)^ 0) ) in = (jM° )f (r > 0)^ 0) ) in = 0, (49) 



Eq. (41), the correlation between '(t) and arbitrary and then the master equation is reduced to 



operator S(t) in Hs system is derived as (p, = c, x) 

(S(t)F(°\t + r)> = - f dt> G„(t + r t')(S(t)s,(t')), 

(45a) 

(FW(t + r)S(t)) = - f dt> G,(t + r t')(s^t')S(t)), 

Jto 

(45b) 



0_ 

dt 



3.k=L,XJ 

+ \phV%\P] 



(50) 



(§(t + T)PW(t)) 



dt' GJt-t')(S(t + T)sJt')) 



■n{5t7(r)[5 M £ difl8 [p(t)] - C diss [3^p(t)}}}, (46a) 



(F^(t)S(t + r))=- / dt' G,(t-t')(s,(t')S(t + r)) 



Tr{SU(T)[C diss [p(t)}^ ~ C diss [p(t)s„}}}, (46b) 



where 



a for p, = c 
b for p = x 



(47) 



Whereas Eqs. (45) are zero in the limit of time-local case 
G M (r) oc S(t) [35], they are in general non-zero in the 
present nonlocal situation. Of course, if r is large enough 
compared to the memory time of the reservoirs, Eqs. (45) 



The steady state obtained from this equation is certainly 
the ground state of closed case p ss = |g)(g|, then the 
decay of polariton system to its original ground state |g) 
is guaranteed by supposing the squeezed and correlated 
reservoirs in Eqs. (48). 

By using this master equation (50), we have calcu- 
lated the dynamics of Hg system. In Fig. 2, supposing 
the vacuum state (no photon and no excitation) at the 
initial time t = 0, the numbers of lower and upper polari- 
tons are plotted as a function of time. In the numerical 
simulation, the density operator p(t) is moved outside 
the time integral, and the memory kernels are also given 
in Eq. (40). While there are non-zero polaritons at the 
initial time, the numbers of polaritons decrease and fi- 
nally go to zero, i.e., the Hg system decays to its ground 
state |g). In the inset of Fig. 2, we also plot the numbers 
of photons and excitations in the early stage. Whereas 
both of them are zero at the initial time t = 0, they 
are oscillated with two periods tt/ujl and tt/lojj (slightly 
modified by the Lamb shifts), but finally they reach to 
(a^a)g = (b'b) g = 0.207 after a long time (not shown in 
the figure). 

Under the Born approximation, the total density op- 
erator /5t t is approximately represented by the product 
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FIG. 3: Sketch of ring-cavity system. Inside the cavity, pho- 
tons interact with matter, but back scattering of photons does 
not occur. The external field is defined in the one dimensional 
system with length L, and the field is continuously connected 
at the boundaries z = ±L/2. 



FIG. 2: Starting from the vacuum state of photons and exci- 
tations, the numbers of polaritons are calculated as a function 
of time by master equation (50) , which are derived by the cor- 
related and squeezed reservoirs as in Eqs. (48). In the inset, 
the numbers of photons and excitations are also plotted in 
the early stage. Parameters: w x = ui c , g = uj c , D = g 2 /ui x , 
\k\ 2 = | 7 | 2 = KTV/27T, and ^ ut " off = ^ ut " off = 10'V. 



of the density operator p of Hg system and the one /3r 
of reservoirs as ptot = P ® PR- If the system-reservoir 
coupling is weak enough for the Born approximation, in 
the ground state |g) of the whole system, the state of Hg 
system is approximately equivalent to the ground state 
g) of the closed case. On the other hand, the free-field 
correlation (48) approximately reflects the reservoir state 
that is obtained by tracing over the Hg variables on the 
ground state |g) as pR ~ Trs{|g)(g|}, which was veri- 
fied in Sec. III. This reservoir state is not the ground 
state of i?R, but it certainly guarantees the decay of Hg 
system to its original ground state |g) as seen in Fig. 2. 
If we suppose the ground state of TJr system, in which 
photonic and excitonic reservoirs are in vacuum (at zero 
temperature), the Hg system does not decay to its ground 
state |g) as seen in Fig. 1. However, the obtained steady 
state is approximately equivalent to the ground state, if 
the system-reservoir coupling is weak enough. Next, let's 
calculate the output from the cavity in the formalism of 
master equation. If the Hg system is in the ground state, 
we cannot detect anything outside the cavity. 

As seen in Fig. 3, we consider a ring-shape cavity em- 
bedding a matter interacting with photons inside the cav- 
ity as discussed in Ref. [35] . We assume that back scatter- 
ing of photons does not occur during the light-matter in- 
teraction, and the clockwise and counter-clockwise fields 
are separated. Concerning the external field, we consider 
a one dimensional system with length L, and the field 
is continuous at the two ends z = ±L/2. Whereas the 
external photonic modes are characterized by wavenum- 
ber kj — 2nj /L for j = 0, ±1, ±2, . . ., the forward 



field j > and backward j < fields can be separated 
into independent subspaccs. Here, we focus on the for- 
ward field j > 0, and its frequency is represented as 
f2j = ckj, where c is the speed of light. The density 
of states (DOS) 

i s n DOS — L/2-kc. This forward field 
couples with the counter-clockwise intracavity photons. 
We define the propagating field in forward direction at 
position z in the external system as 

F c fwd (M) = J2 K ma m (t)c ik - z . (51) 

/ n 

From the equation of motion of a m (t), this field is rewrit- 
ten as 

m 

+ f dt' G c (t-t' - z/c)a(t'). (52) 

Jta 

As discussed in Ref. [35], by choosing a position of ob- 
servation z > 0, we define the output field as 

P^(t)=F^ d (z Q ,t + z /c) 

rt + z /c 

= H 0) (t) + / dt' G c (t - t')a{t'). (53) 

Here, Fc°\t) is the free field appearing in the Langevin 
equation (41) and also in the master equation (31). The 
second term is the contribution from the cavity. Whereas 
this term includes the information of cavity photons at 
time t + zq/c > t, the causality is not violated, because 
the output field F° ut (t) is actually the propagating field 
at position zq and at time t + zq / c. In the time- local limit 
G c (t) = r c S(r), Eq. (53) is correctly reduced to the well- 
known input-output relation [35, 36, 38]. Further, in the 
limit of zq — > oo and t — > — oo, Eq. (53) is reduced to 
the input-output relation (60) in the time-nonlocal case, 
which will be derived in Sec. V. 
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When we evaluate the output measured by photon de- 
tectors, the expectation values should be normal-ordered 
and time-ordered (expressed as (: ... :)) in terms of po- 
lariton operators (not of photon and excitation). The 
correlation between cavity photons and the free field of 
photonic reservoir can be evaluated by Eqs. (45) and 
(46), and the self-correlation of the free field is also given 
by Eqs. (48). The detail of the calculation is shown in 
App. C. When we suppose that the H$ system is in the 
ground state |g) by considering the reservoir correlation 
(48) in the master equation, we have numerically checked 
that the emission spectrum (: F° nt (w) F° ut ^ :} and phase- 
sensitive correlation (: F° ut (uj)F° ut :) are approximately 
zero. The deviation is due to the approximation that we 
used in the numerical calculation (density operator p(t) 
is moved outside the integral) , and it is not caused by the 
supposed correlation, Eqs. (48). On the other hand, if we 
suppose the vacuum photonic and excitonic reservoirs, we 
cannot find a policy which guarantees no photon detec- 
tion, although the vacuum output is obtained for vacuum 
input in the input-output formalism [2] . This is because 
of the perturbation treatment in the formalism of master 
equation as we will discuss in Sec. VI. 

In this way, when we suppose the squeezed and cor- 
related reservoir fields as in Eqs. (48), we have success- 



fully obtained the natural result: the Hg system decays 
to its ground state |g), and the photon emission is not 
detectable if the system is in the ground state. Further- 
more, it is also consistent to the analysis of Fano-type di- 
agonalization (there are virtual photons and excitations 
in the reservoirs, and photonic and excitonic reservoirs 
are correlated with each other and also squeezed). 



V. INPUT-OUTPUT APPROACH 

Another approach for describing the dissipation and 
emission of photons is the formalism of Langevin equa- 
tions with input-output relation. As discussed in Ref. [2] , 
the Langevin equations of cavity photons and excitations 
arc derived in frequency domain as 









S(w) 


= i 




o(-w)t 


F^ n (-w)t 


6(-w)t 




J>(-a;)t. 



[M(w) -wl] 



Here, the coefficient matrix is written as 



(54) 



M(w) 



w c + 2D — iG c (uj) + ig 

-ig oj x -iG x (uj) 
-2D -ig 







2D 
-i.9 

-u c -2D- iG c (-w) 
-i.9 



-1.9 


-cj x - iG x {-uj)* + . 



(55) 



and the memory kernels G c _ x (t) arc Fourier-transformed 
for positive time as 



G c , x (u)) + = / dr e lUT G CtX (T). (56) 
Jo 



The Langevin (fluctuation) operators are expressed as 



Their Fourier transforms are derived as 



(*-*<>) 



(57a) 

£j» = E^™(*o)e- ifi - (t -* o) = $>w3£e- ifi ™*, 

(57b) 



where to — > — oo is the switch-on time of system-reservoir 
interaction, and a™ and /3™ are the input operators. 



dt e iu]t F^{t) = 2Tr6{aj) K (u)a in (uj), 



(58a) 



dt e iu]t F x n (t) = 27r6»(w) 7 (u;)/3 in (c l ;). 



(58b) 



Here, the reservoir states are rewritten in continuous 
form as in Eqs. (Al). These fields F™ x (uj) are interpreted 
as the input fields, and they cannot be defined for neg- 
ative frequency uj < 0, because the reserver states are 
distributed only for positive frequencies ff^ x > 0. 

According to the input-output formalism [2], the out- 
put photonic field (photonic reservoir field at time t\ — > 
oo) is represented as 



d out (w > 0) = d in H + k(w)*o(w). 



(59) 



As discussed by Ciuti and Carusotto [2], we get the 
vacuum output for vacuum input. However, the Hs 
system is actually excited by the vacuum reservoirs as 
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(Pj(u)pk) 7^ and (pj(uj)pk) ^ 0, which can be easily 
verified from the Langevin equations (54). In the master- 
equation formalism discussed in the previous section, the 
Hg system is also excited, but the vacuum output is not 
obtained for the vacuum input. Then, there is a discrep- 
ancy between the two approaches at least under the Born 
approximation. Instead, in the input-output formalism, 
we also suppose the squeezed and correlated reservoirs 
discussed in Sees. Ill and IV. 

According to the standard input-output formalism, the 
output photonic field is represented as 

/oo 
dt' G c {t - t')a{t'), (60a) 
-OO 

F c out M = ^ n H + G c Ha( W ). (60b) 



This expression does not violate the causality (F° ut (t) 
can be affected by a(t' > t)) as discussed in Sec. IV. 
From this input-output relation and the Langevin equa- 
tions, the output photonic field is eventually represented 
by the input fields F™ x (uj). For discussing the output 
from the cavity, we have to suppose the correlation of 
input operators {F™ x (uj)}. Here, we consider that the 

Hs system is in the ground state, and the correlation of 
input operators are also supposed as shown in Eq. (48): 



(61) 



(i>H + i^ T ) R = 



/ G c (cj) + (aa t ) g G c (w)+(a& t ) g 
G c {-u)* + {tftf) & G c (-Gj)* + {tftt) g 

\G x (-w)5.(Stot>g G x (-uy + (W) t 



G c (uj) + (aa) g G c (u) + (ab) g \ 

G x {uj) + (ba) g G x (uj) + (bb) g 

G c (-w);<a+a) g G c (-w);(a+6) g 

G x (-u)* + (tfa) g G a (-u)* + ftb)J 



/G c H;(aat) ( 
G c (u))* + (tftf) 



G x {Lo)* + {ab^) i 

G x (w)* + {aW) 
G x (wy + {bW) : 



G c (-uj) + {aa) g 
G c {~u) + {ba) g 
G c (— Lu) + (a^a) g 
G c (-a;) + (Sta) g 



G x (-uj) + (ab) g \ 
G x (-iu) + (bb) g 
G x (-w) + (a t &) g 

G x (-Lj) + (ttb)J 



where Pg(u) = [P?(u), F™(uj), F^-lj)^ ^(-w/f. 
Precisely speaking the expectation values such as (a^a) g 
should be slightly modified depending on uj as discussed 
in App. D. Assuming this input correlation, the system is 
certainly in the ground state (pj(uj)pk) = (pj(oj)pk) = 
(j, k = L,U), and we also obtain no photon detection 
(: F° ut {uj)F° u ^ :} = and (: F c out (w)F c out :) = by 
normal- and time-ordering the operators in the polariton 
base. The detailed calculation is shown in App. D. 

In this way, when we suppose the squeezed and cor- 
related reservoirs represented in Eqs. (48) and (61), the 
system certainly decays to its ground state |g) and no 
photon is detected outside the cavity in both formalisms 
of master equation and input-output relation. In con- 
trast, when we suppose the vacuum reservoirs, different 
results are obtained in the two formalisms. 



VI. DISCUSSION 

As discussed in the previous sections, when we con- 
sider the squeezed and correlated reservoirs instead of 
the vacuum ones, both master-equation and input-output 
formalisms certainly guarantee the decay of the Hg sys- 



tem to its ground state and show no photon detection, 
even though we do not use the RWA on system-reservoir 
coupling. The supposed reservoir state is approximately 
equivalent to the one realized in the ground state of the 
whole system: p R ~ Tr#{|g)(g|}. On the other hand, if 
we suppose the vacuum reservoirs in excitation-photon 
base, in the absence of the RWA on system-reservoir 
coupling, the Hg system is excited by the coupling with 
the reservoirs as seen in Fig. 1. We have to determine 
the reservoir state supposed in the master equation and 
input-output formalism, according to the situation how 
the system and reservoirs start to couple. If we initially 
prepare the vacuum reservoirs and switch on the system- 
reservoir coupling, the reservoirs approximately remain 
in the vacuum state even after the switch-on, and the 
system does not decay to its ground state but to a steady 
state excited by the vacuum reservoirs as seen in Fig. 1. 
This is because the Hg system is excited when virtual 
photons escape to the vacuum reservoirs. In order to 
avoid it, the -Hs system and the reservoirs should be bal- 
anced as realized in the ground state |g) of the whole 
system, and we should suppose /Sr ~ Trjj{|g)(g|} in such 
situation. If we consider that the system and reservoirs 
are already coupled and the whole system is in the ground 
state, when we excite the H$ system to an excited state, 
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the system certainly decays to its ground state as seen in 
Fig. 2. 

If the reservoirs are quite large and the whole system 
cannot be in a steady state, we should suppose the former 
situation. When the temperatures of the reservoirs are 
low enough and the vacuum input from the reservoirs to 
the system is supposed, the Hg system in principle does 
not decay to its ground state. Although the vacuum out- 
put is obtained according to the input-output formalism 
[2], it is not by the master equation. The energy is con- 
served in the input-output formalism, but it seems not 
in the master equation. This is because the dynamics 
of focusing system and the output are discussed in the 
sense of perturbation theory in the formalism of master 
equation. In this way, when we suppose the vacuum in- 
put, we should pay attention to the difference of the two 
formalisms (at least under the Born approximation). In 
order to avoid this discrepancy, we should use the RWA 
on system-reservoir coupling, although the quantum fluc- 
tuation of reservoirs is diminished in such treatment. 

On the other hand, if we can define relatively small 
reservoirs which weakly couple with a large external sys- 
tem with low enough temperature, the small reservoirs 
and the Hs system can decay to the ground state |g) of 
the coupled system. In such situation, we can suppose 
Pr ~ Trfl{|g)(g|}, and it guarantees the decay of Hs 
system to its ground state |g) and gives no photon de- 
tection in the small reservoir as discussed in the previous 
sections. This result is obtained in both formalisms of 
master equation and of input-output relation in contrast 
to supposing the vacuum reservoirs in excitation-photon 
base. 

As discussed in Ref. [24], by performing the RWA on 
system-reservoir coupling, we can simply suppose the 
vacuum reservoirs in excitation-photon base, and the 
master equation is reduced to the standard Lindblad 
form. The simplified master equation is derived as fol- 
lows. Whereas the system-reservoir coupling is originally 
represented as Eq. (9), here we perform the pre-trace 
RWA [41] as 

H S . R ~ ih E {WjF^pj ~ Wjp]P c + x*F^pj - Xjp)P x ), 

j=L,U 

. (64) 

where the counter-rotating terms are neglected in the 
polariton base not in the excitation-photon base. Then, 
when we suppose the vacuum reservoirs in the excitation- 
photon base as in Eqs. (39), the master equation is de- 
rived under the Born approximation as 

+ E / dt ' {Gc{t-t')w J wl + G x {t-t')x j x* k } 



j,k=L,U^ l ° 



[u s (t-t')\p k p(t%p\ 



H.c. 



(65) 



Further, by neglecting the fast oscillating terms 
-i(,u k - Uj )t for j ^ k (called the post-trace RWA 



[41]), wc finally get the simplified master equation under 
the Markov approximation as 

d - 1 in -i 



3=L,U 



J 1 

(66) 



where the memory kernels are approximated as G M (i) = 
r^Slt) for simplicity (there remains the Lamb-shift terms 
in general [24]). If the system-reservoir coupling is ex- 
pressed in the Hcrmitian form as ih(d ± aJ)(FJ =F F c ) 
[ih(b±tf)(F%=fF x )], the above master equation is simplify 
rewritten by replacing wj [xj] by Wj =Fj/j [xj^Zj]. Even in 
such case, the simplified master equation is represented 
in the Lindblad form. From Eq. (64), the input-output 
relation is obtained as 



F° 



F" 



(67) 



Since the above master equation is reduced to the stan- 
dard form owing to the pre-trace and post-trace RWAs, 
we consider that the correlation of input operator i 7 "" 1 

is equivalent to that of F^ supposed in the master 
equation: (F™(t)F^) R = G c (t) and (F^(t)F™) K = 
(Fj n (t)# c in ) R = 0. Then, the correlation of the output 
can be calculated as discussed in Ref. [28]. 

However, in this approach, the photonic and excitonic 
reservoirs are supposed to be in the vacuum state un- 
der the RWA on system-reservoir coupling, although the 
polariton system does not decay to its ground state |g) 
in general without the RWA. In other words, the quan- 
tum statistics of reservoirs fields are diminished by the 
RWA, although the reservoirs are originally squeezed and 
correlated. In contrast, in the present paper, the master 
equation and input-output formalism arc discussed based 
on the squeezed and correlated reservoirs. The master 
equation certainly guarantees the decay of H$ system to 
its ground state, and in both formalisms no photon is 
detected when the system is in the ground state. Un- 
der the Markov approximation the master equation (50) 
is reduced to 

q 1 - F- 



+ E 

j,k=L,U 



j,k=L,U 

^7T- (PjPPk - PkPjP) + H.c. \ . (68) 



PkPjG 



The coefficients and Kj^ can be calculated from the 
supposed free field-correlation in Eqs. (48). This does 
not have the Lindblad form, but certainly guarantees the 
decay to the ground state |g) as seen in Fig. 2. 

In the standard theory [36-38] and also in the dis- 
cussion of Refs. [24, 28], the master equation and the 
input-output relation are sometimes used together and 
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the correlation of input F™ is supposed to be equal to 
that of free field given in the master equation. How- 
ever, the formalism of master equation is discussed in 
the sense of perturbation theory. Since the reservoirs are 
large enough compared to the H§ system, the input cor- 
relation is not strongly modified and constantly given in 
the master equation. On the other hand, the output is a 
perturbation of the reservoirs as a result of the system- 
reservoir coupling. The correlation of Fc°^ can be in gen- 
eral different on input and output sides. Actually, when 
we suppose the vacuum reservoirs in excitation-photon 
base, the self-correlation of Fc ^ on output side is not 
in vacuum, which is calculated by Eq. (48). In order 
to get the same correlation for input and output sides, 
we have to consider the squeezed and correlated input 
/5r ~ Ti'j}{|g)(g|}. If we want to reduce this complicated 
formalism into the standard one, we have to perform the 
RWA on system-reservoir coupling [24, 28]. 

If we already know that the free field F^ does not con- 
tribute to the observables, we can simply use the RWA 
on system-reservoir coupling [24, 28]. For example, the 
second-order correlation functions under a resonant exci- 
tation can be calculated as discussed in Rcf. [28]. How- 
ever, when we discuss squeezing of the emission, the in- 
terference between free field Fc°^ and cavity contribution 
is important, and the quantum fluctuation of Fc°^ should 
not be destroyed by the RWA on system-reservoir cou- 
pling. If the cavity system has an optical nonlinearity or 
embeds ensemble of atoms, we have to treat the Langevin 
equations pcrturbatively or the master equation might be 
appropriate to treat such systems. When we discuss the 
emission (or lasing) from such complex systems under in- 
coherent excitation, it is difficult to evaluate the validity 
of the RWA on system-reservoir coupling, and we should 
suppose the squeezed and correlated reservoirs realized in 
the ground state of the whole system. This kind of ap- 
proach should give us natural results in the calculation 
of dissipation and detection of output. 



VII. SUMMARY 

We have derived the master equation, Langevin equa- 
tions, and input-output relation for dissipative polariton 
system in the ultrastrong light-matter coupling regime. 
The correlation of reservoir free fields are required for 
calculating not only the dynamics of the system but also 
the photon emission from the polariton system. When 
we suppose the vacuum reservoirs, the polariton system 
is excited in general. Although the vacuum output is ob- 
tained for the vacuum input in the input-output formal- 
ism, it is not obtained in the master-equation approach 
under the Born approximation. In order to avoid this 
discrepancy, we have to perform the RWA on system- 
reservoir coupling, although it diminishes the quantum 
statistics of the reservoirs. In order to describe the dis- 
sipation in the ultrastrong coupling regime without the 



RWA on system-reservoir coupling, we have considered 
the correlation functions of the photonic and excitonic 
free fields that are squeezed and correlated with each 
other and realized in the ground state |g) of the whole 
system: /5r ~ Trg{|g)(g|}. In the formalism of master 
equation, the supposed correlation certainly guarantees 
the decay of the polariton system to its original ground 
state |g). In the ground state, we have also verified no 
photon detection as the output from the cavity. Even in 
the formalism of Langevin equations and input-output re- 
lation, we also get no photon detection by considering the 
squeezed and correlated reservoirs. This reservoir state 
is also consistent to the analysis of the ground state of 
the whole system by the Fano-type diagonalization tech- 
nique. At least when the polariton system is dissipative 
and is in the ground state, the three approaches, master 
equation, input-output formalism, and Fano-type diago- 
nalization give the same result, in contrast to supposing 
the vacuum reservoirs. The case in the presence of exci- 
tation to the system will be discussed in the future. 
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Appendix A: Diagonalization of photonic and 
excitonic parts 

In order to diagonalize the whole Hamiltonian H = 
Hs + ffs-R + -ffRj first of all, we diagonalize the photonic 
part, Eq. (12). Here, we rewrite the reservoir fields from 
discrete to continuous form as 

(ffa), (Ala) 

p m Pm/JnZosinz), (Alb) 

K m -)• ^J/y^og^), (Ale) 

7m 7(ft£)/^n£os( fi m)- ( Ald ) 

where nQ s( w ) an d ^dos^) are densities of states 
of photonic and excitonic reservoirs, respectively. 
The new reservoir operators satisfy [a(ui), d T (a/)] = 
[(3(uj), p'{u)')\ = 6(u> — ui'). The photonic Hamiltonian 
is rewritten as 

/>oo 

Hph = kuj c a)a + / duj ftwa^ (uj)a{<jj) 
Ja 

/>oo 

+ ih dw [«(o;)*d t (w)a - a t d(w)K(a;)] . (A2) 
Jo 
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This kind of Hamiltonian can be diagonalizcd by the 
Fano-type technique [31-33] by introducing an operator 
for eigen frequency u> as 

/>oo 

A(u) = u c (uj)a + / du/ v c (uj 7 uj')a(uj'). (A3) 
Jo 

Once this operator satisfies Eq. (13), we can diagonalize 
the photonic Hamiltonian as in Eq. (12). Further, A(uj) 
should be normalized as 



AH,AV) 

poo 

= u c (uj)u c (uj')* + / duj" v c (uj,uj")v c (uj',uj")* 
Jo 

= 6(uj-uj'). (A4) 



The coefficient functions u c (lu) and v c (oj, uj') are deter- 
mined as follows. From Eq. (13), we get 



where 



uju c (uj) = uj c u c (uj) + i / du/ v c (uj, uj')k(uj'), (A5a) 
Jo 

uiv c (u!,lu') = uj'v c (uj,uj') — ik(uj')u c (uj) . (A5b) 
From the second equation, v c (uj,uj') is expressed as 
v c (u,u') 

(A6a) 



-ik(uj ) — 



-i/c(u/) j P + ib(u;)5(uj - w') \ uAuj), (A6b) 

[ UJ — uj' J 



-Ik(uj') 



-uj 1 - i0+ 



+ ["0( w ) — m]5{uj — uj') f u c (uj). 

(A6c) 



where P means the principal value integral and function 
ij>(u)) is introduced for the following calculation. The ex- 
pression of is determined by substituting the second 
or third equation into Eq. (A5a) as 



i>{uj) 
"4>{uj) — \n 



1 



1 

WW 



duj 



A|2 



(A7a) 



w-u)' - i0+ 



(A7b) 



C(w) = 1 - - r da/ J^ 01 ' . (A10) 
v ; uj c J uj' -cj + iO+ v ; 



Using the diagonalized operator A(uj), the original ones 
are represented as 



= / duj u c {uj)* A{uj), 
'0 



(Alia) 



a(w)= duj 1 v c (uj',uj)*A(uj'). (Allb) 
Jo 

In the same manner, we can also diagonalize the exci- 
tonic Hamiltonian as in Eq. (15). The eigen operator is 
represented as 

/>oo 

B(uj) = u x {uj)b+ duj' v x (uj,uj')(3(uj'). (A12) 
Jo 

The coefficient functions are determined in the same 
manner by replacing w c and k(oj) with uj x and "/(uj), re- 
spectively. The excitations and excitonic reservoir field 
are represented as 



du u x (uj)* B(uj), 



(Al3a) 



(3(uj)= duj' v x (uj' ' ,uj)*B(uj'). (A13b) 
Jo 



Appendix B: Correlation of free field 

By using the technique in Ref. [35], here we calculate 
the correlation between the free field Fc° x (t) and sys- 
tem operators a and b. Further, we also derive the self- 
correlation of Fc° x (t) . 

First, let's calculate (S(t)Fc (t + r)) for arbitrary sys- 
tem operator S and t > 0. From the Langevin equation 
(41), the free field is represented as 



*f>(t) 



d , 1 



a(t),H £ 



then we get 



d*' G c {t-t')a(t'), 
(Bl) 



On the other hand, the expression of u c (uj) is determined 
by the normalization condition, Eq. (14). The commuta- 
tor is derived as 



A(oj),A^(u}') = u c (uj)u*(uj')[^(uj)-m}[^(uj)+m}8(uj-uj'). 

(A8) 

Then, we get 



{S(t)F^(t + r)) 

= ~0(t)a(t + r)> + ^(S(t)[a, H s }(t + r)> 



t+T 



di' G c (t + T-t')(S(t)a(t')). 



(B2) 



la 



u c (uj) 



\n(uj)f 



iJj(uj) — m uj — uj c ^(uj) ' 



(A9) 



The correlation functions of system operators appearing 
on the right hand side can be calculated by the master 
equation (31). By using the quantum regression theorem 
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(44), the first term of Eq. (B2) is rewritten as 



and we get Eq. (45b). 



= Tr{a£[U(r)[p(t)S}}} 

= ±(S(t)[a,H s ](t + T)) +Tr{al diss [U(r)[p(t)S]]}, 

(B3) 



and its last term is also written as 

j-t + T 

Tr{a£ diss [U(T)[p(t)S}]} 



dt' (S(t)a(t'))G{t-t r ) 
(B4) 



The next is (S{t + T)F^ 0) (t)) for r > 0. From Eq. (Bl), 
we get 



= -(S(t + r)^a(t)> + -<£(* + r)[a,H s ](t)) 



dt' G c (t - t')(S(t + T)a{t')). (B7) 



to 



Substituting these two equations into Eq. (B2), we get 
Eq. (45a). 

Next, let's consider (PP(t + r)S(t)) for r > 0. From 
Eq. (Bl), we get 

(F^(t + r)S(t)) 

= - J^(* + r)S(t)) + r\{[a(t + r),H s }S(t)) 



t+T 



dt' G c (t + T-t')(a{t')S(t)). 



(B5) 



to 



As shown in Rcf. [35], the first term is represented as 



(S(t + r)-a(t)) 

= l(S(t + r)[a,H s }(t)} 
in 

+ Tr{SU(T)[a£ diss [p(t)} - AlissM*)]]}, (B8) 



In the same manner as the above calculation, the first 
term is rewritten as 

jL{a(t + r)S(t)) 



t+T 



and then we get Eqs. (46). The second equation is also 
derived in the same manner, and similar expressions are 



obtained also for F 



(0) 



-([a,H s ](t + r)S(t))- jf dt' G c (t + r - t'){a(t')S(t)}, oj 



Since the free field Fc°\t) is expressed as in Eq. (Bl), 
(B6) the self-correlation is represented as 



<Fj >(i)Fj°>t( T )> = JLJL( 6 ( t )at (T)) _ ±l{a(t)[tf,H B ](r)) + f dt' G* c (r - t')^-(a(t)tf(t')) 



m dt 



to 



i__a_ 

ih dr 



[a,H s )(t)aHr)) + —r ~([a, H a ](t)[tf , H 8 ](t)) 



(ih) 



dt 



dt' G* c (r-t')([a,H s }(t)tf(t')) 



to 



+ f dt" G c (t-t")^-(a(t")a^T)) 

J to dT --ro 



lS Jtr 



dt" G c (t-t"){a(t")[a\H s }(T)) 



dt" dt' G c (t-t")G*{T-t')(a{t")a\t')). (B9) 

J t 



In the same manner as discussed above, the first term is rewritten for r > t as 



||-(a(t)at(r)) = ^Tr{a^U(r - t)[p(t)a}} 



I|-(a(t)[at,F s ](r))-| / 



ih dt 
1 d 



^{a{t)[a\H s ](r))- f dt' ?-(a(t)tf (t'))G* c (r - 1>) + G* c (r - t)(a(t)a!(t)}. (BIO) 
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Rewriting 4th and 7th terms in Eq. (B9), the self-correlation is reduced to 
(Fi°)(t)Fi°)t(r)) = G c *(r - t)(a(t)aHt)) + J* dt' G*(r - t')^{a{t)a\t')) 

&t' G* c {T-t'){[a,H s ]{t)a){t'))+ ( dt" f dt' G c (t - t")G*(r - t'){a{t")a) \t')) . (Bll) 

WJto J t J t 

Rewriting the second term, we finally get 

(F^(t)F^(r))=G* c (r-t)(a(t)a t (t))- f dt' f dt" G c (t t")G*(r - t'){a{t")a ] '(*')> 
+ f dt" f dt' G c {t-t")G* c {r -t')(a(t")a){t')) 

J t J t 

= G*(r-t)(a(t)a\t)). (B12) 
In the same manner, we finally get Eqs. (48). 

I 



Appendix C: Calculation of observables in 
master-equation approach 

When we detect photons emitted from the cavity, the 
observables by photon detectors should be calculated by 
normal- and time-ordering the photon operators [35, 36, 
38] . In the present case, the ordering should be performed 
in the polariton basis, which really represents the eigen 
states of the system. Here, we divide the photon operator 
a into the lowering parts fij, and raising part as 

h = w* L p L + w^Pu, (Cla) 

Of = -VlPl - VuPu- (Clb) 

Since the system-reservoir coupling is expressed as in 
Eq. (9), the photonic free field is also divided as 

F<®= W IFM +W * u Pj J °\ (C2a) 
F^=-y L F^- yu P^, (C2b) 
and then the output field (53) is rewritten as 

pout = pout + pout ^ (C3) 



ft+zo/c 

FTV = F S\t) + / & G c (t ty H {t'), (C4a) 

Jto 

pt+zo/c 

F°f(t) = F$>(t) + / dt' G c (t t')a t (t'). (C4b) 

Jto 

The equation of motion of aj,-)- is derived as 

J^t(*) = ^&t.#s](*) 

- f dt' G c (t-t')a it (t')-F%l(t). (C5) 



Then, the correlation between fij.f and F^ is also derived 
in the same form as Eqs. (45) and (46), and the self 
correlation of F^l also has the same form as Eq. (48). 

The output field measured by photon detectors outside 
the cavity is calculated by normal- and time-ordering the 
operators. The emission spectrum (number of photons) 
in a steady state is expressed as 

(: Pr\uj)F^ :) 

= (:FW(uj)FW :) + G c (u)(:&(uj)FW :> 

+ G c H*(:F c (0) Hat :) + \G c (uj)\ 2 (: a(w)a t :), (C6) 

where 

(: Fj°>( W )Fj°)t :) 

= G c (cj){(a|%) + (a t a{)} + G c (w) + {(a|a t ) + (a^aj)} 
+ G c (u)X{(a t a\) + (a\a i .}}, (C7a) 

(: a(u>)FW :) 

= (f$\ u )%) + (f^iwt&r) + (« t (-)-^ 0)t ) 

+ <a,H + ^ T 0)t ) + (F^i^la,) + (a t ^)f c ( t 0)t ), 

(C7b) 

<:F c (°)(a;)at :) 

+ (F^(uj) + a\) + (a t (a;)t#S ) ) + (F c ( t 0) ^)a|), 

(C7c) 

(: a(uj)a^ :) 

= (o 4 .(w) t o 4 .> + (a4.(w)Lat) + (a t (a;)_a{) 
+ (a^(u;) + a{) + <a t M+«i) + <«t(w)4)- ( C7d ) 
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On the other hand, the phase-sensitive correlation is ex- 
pressed as 

(: F° ut (w)F° ut :} 

= (: PM^PM :)+G c H<:aHFj°> 
+ G c (-cu)(: F^(tu)a :) + G c (lu)G c (- 



(:^ (0) M^ (0) :) 

= {Gc(u)+ + G c (-uj)+}{2(a t a i ) + (a t a t ) + (aj.%)}, 

(C9a) 

(: o( W )Pf) :) 

= (a^+F^) + (F^i-uUa,) + (F§\( U )) 

+ (a t (a;)#ip) + (F c ( t 0) (- W )_a t ) + (a t (w).F^}, 

(C9b) 

<:Fj°)Ha:> 

= (F c ( , 0) ( W ) + a i ) + (a,(- W ) + F< i 0) ) + <M^%)> 

+ (F^(w)aj.) + (a t (- W )_Fi t 0) ) + (F< t 0) (a,)_o t >, 

(C9c) 

(: a(w)d :} 

= (%H+ a 4-) + (a-li-^+ai) + (a t (-u;)%) 

+ (a^(u)a±) + (af(— a;)- Of) + (a t (a;)_a^). (C9d) 

We have numerically checked that the correlation func- 
tions (: #° ut (w)F° ut l :) and (: F° ut {w)F° ut :} are ap- 
proximately zero if the Hg system is in the ground state 
p ss = |g)(g|. The small deviation comes from the ap- 
proximation in which the density operator pit) is moved 
outside the time integral in the master equation (31). 



Appendix D: Calculation of observables in 
input-output formalism 

Let's calculate the emission and squeezing of photonic 
output from the cavity in the input-output formalism. 
The coefficient matrix (55) of the Langevin equations is 
diagonalized as 



M(oj) = V(w)D 



u>u(cj) 
-Q L (-u)* 
-Qxj{—u})* 



v(«)- 



(Dl) 



where D[- • • ] represents an diagonal matrix with elements 
• • ■ . Due to the coupling with reservoirs, the eigen val- 
ues {wi.[/(w)} depend on frequency uj, and are modi- 
fied from the original eigen frequencies {u)l,u} derived 



from Eq. (4) . The modification depends on the strengths 
Kj and 7j of system-reservoir coupling. Wc redefine 
Langevin (fluctuation) operators in the polariton basis 
as 



) 








r # o) m " 




a(uj)a :), 
(C8) 






= v( w )- 1 


^i o) (^) 
^ (o) (-^)t 

F^(-u>V 


■ (D2) 



Because of the modification of the coefficients, these op- 
erators are in general different from the Fourier transform 
of free field of the reservoir field Pl,u, Eq. (10), in the po- 
lariton basis. However, if the system-reservoir coupling 
is weak enough compared to the characteristic frequency 
of the polariton system, the redefined operators are ap- 
proximately equal to the Fourier transform of Eq. (36). 
At the same time, the Born approximation used in the 
master equation is also valid. The photon and excitation 
operators are then represented in the frequency domain 
as 



a(oj) 
b(uj) 
o(-w)t 

5(-w)t 



where 



L(w) = [M(w) -wl] _1 V(a;). 



(D3) 



(D4) 



Substituting Eqs. (D2) and (D3) into it, the photonic 
output operator is represented as 



Fr M = T L (c)F { L a) (c) 

*j L °\-«,)t + S u (u>)P<> 

where 



■AO), 



+ S L (w)Pi°\-w)* + Su^P^i-u^, (D5) 



Tl(lo) = V n (cj) + G c (w)L u (w), (D6a) 

Tu(lu) = V 12 (lu) + G c (lu)L 12 (lu), (D6b) 

Sl(lo) = V 13 (lu) + G c (lu)L 13 (lu), (D6c) 

Su(lu) = V u (lu) + G c (u))L 14 (uj). (D6d) 

As discussed in Sec. V, for describing the dissipation 
of Hs system, we consider the correlation of free fields as 
in Eq. (61). Precisely speaking, in order to guarantee no 
photon detection, the expectation values in the ground 
state should be replaced by 



/ (aa f )g 
(Sot) g 

<atat) g 
V<6t a t> g 



{ab% 

(atSt )t 



(aa) g 
(So) g 

(Sta) g 



(ab) B 

(SS) g 



) 


<l 










-».vh 




i 



















v!! 








0/ 



V( 



(D7) 

because of the modification of coefficients V(w). If we 
detect the output photons outside the cavity, the detec- 
tion process should be dissipative for the whole system, 
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and virtual photons should not be counted. Then, when 
we calculate emission spectrum (: F° ut (w) F° ut t :) and 
phase-sensitive correlation (: F° ut (cu)F° ut :). the fluctua- 
tion (Langevin) operators {f£ ^(w)} should be normal- 



and time-ordered obeying the theory of measurement 
[35, 36, 38]. Then, the correlation functions of output 
photonic field (D5) are represented as 



{: F^i^F^ :) 



j,k=L,U 



-S» [(F^i-^lF^) + (Ffi^lF^)] T k {uf + 5 J -(^)(Fj 0) (-a;)+F fe (0) )5 fc (a;)*} , (D8) 



(: C'(w)i ? r :> 



j,k=L,U 



- u 

+S j (u>)(Fj°\-^Fl; 0) )T k (-cj) + S j (uj) [<ff (- W )!ff> f ) + (if W-^f f }] S h (-u>)} . (D9) 



If we consider the fluctuation correlation as shown in are completely zero. 
Eq. (61), we can numerically verify that both of them 
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